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Abstract 

We study the simultaneous influence of boundary conditions and external fields on quantum fluctuations 
\q ' by considering vacuum zero-point energies for quantum fields in the presence of a magnetic fluxon confined 

by a bag, circular and spherical for bosons and circular for fermions. The Casimir effect is calculated in a 
generalized cut-off regularization after applying zeta-function techniques to eigenmode sums and using recent 
techniques about Bessel zeta functions at negative arguments. 

1 Introduction 

• rH 

x 

Aharonov-Bohm H] settings may be regarded as one of the possible ways in which an external field modifies some 
observables of a given quantum system. In this specific phenomenon, the presence of an infinitely thin tube of 
magnetic flux alters the energy spectrum and brings about a modification of the vacuum energy, giving rise to 
a form of Casimir effect. Initially, Aharonov-Bohm fields acted on free particlesQ The purpose of the present 
paper is to study the influence of the same type of fields on systems which are already constrained by boundary 
conditions (b.c), working out the combined net effect of both on what would otherwise be a free system. The 
relevance of the Aharonov-Bohm scenario to some cosmic strings models including particles in the gravitational 
field of a spinning source is discussed in [| . A similar mathematical procedure is also applied to the description 
of Dirac fermions on black-hole backgrounds as shown in Q . 

The quantum mechanical problem of a scalar particle inside a circular Aharonov-Bohm quantum billiard 
H"H (°f radius a) bears a great resemblance, from the mathematical point of view, to the ones we set out to 
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consider. In our case we have a classical magnetic fluxon which is coupled to a quantum field. We take this 
object to be an idealization of a vortex with a radially symmetric distribution of magnetic field in the limit where 
its characteristic thickness is vanishingly small. With this model in mind one has the physical basis to fix the 
boundary conditions at the origin to be imposed on the modes of the matter field, (a detailed account of this 
kind of analysis is given in pj). We start with a complex, Klein-Gordon, massless field (the additional analytic 
effort which takes the treatment of a massive field by zeta function techniques is explained in We call </> 

the space-dependent part of the eigenmodes, which satisfies the equation (in units such that Ti = c = 1) 



-iV - el) cj) = uj 2 (j), (1.1) 



where the vector potential A is given by 

01 „ e$ ,„ „. 

eA(r) = -e v , ot=—. 1.2 

a is called reduced flux, being $ the flux of the magnetic field. Since a billiard is a domain with perfectly reflecting 
walls, and we imagine an infinitely thin solenoid at the origin — reduced, in D = 2, to an unreachable point — 
the b.c. are <f> = at r = and r = a. 

Zero-point energies emerge from mode-sums of the type — ^^uj n , and give rise to the Casimir effect JT^-Jl^]. 

n 

Since the summation extends over all the w n 's in the set of eigenmodes, such quantities usually diverge and need 
some regularization to make sense of them. To this end, we introduce spectral zeta functions as mere auxiliary 
tools, which will be denoted by 

n n ^ ' 

fi is an arbitrary scale with mass dimensions, used to work with dimensionless objects. This is a regularization 
of analytical nature (see also p^|); other examples in this same category are the techniques in refs. [[l5| and refs. 
|fl6|| , [[[7). When we discuss the results, we shall comment on their physical significance from the perspective 
of cut-off regularization. In this sense, our standpoint in the present case is that the zeta function is a purely 
mathematical object which affords a convenient method for the calculation of observables inasmuch as it may 
be connected with other, more physical, regularizations. Let us assume that we are using a general cut-off 
regularization for the vacuum energy which is given by 



-&reg 

where g is a well-behaved function which satisfies asymptotic expansions near the origin and at infinity of the 
following kind: 

00 1 oc 

g(t) ~ 1 + £ a k t k (t -> 0) g(t) ~ ^ £ b k r k (t -> 00). (1.5) 

k=l k=Q 

If we restrict this analysis for the sake of simplicity to a (2 + 1) — D case, then g should be such that M s > 3. This 
conditions guarantee that the Mellin transform of g has no poles in the strip < SJfz < M s . The (^-function, 



defined in (1.2), has its rightmost pole at s = 2. In fact, the requirement which we demand on the asymptotic 



behaviour of g might have been stated with greater generality. Essentially, g should go to 1 for small values of its 
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argument, and should go to zero for large values fast enough so that the following reasoning makes sense. Now we 



explain presently how the connection between the Cx-function (1.2) and expression ( |1.4| ) may be accomplished 
(for other discussions on how the results from different regulators may be related, see 
The Parseval formula for Mellin transforms allows us to write 



£re S (A) = i-^/ dzk 1+z M[g,z + l]C M (z), (1.6) 



r+zoo 



/ r — zoo 

where r is any real number such that 2 < r < M s — 1. The structure of the poles of the integrand in expression 



(1.6) is closely related to the asymptotic expansion of E reg (A) for large values of the cut-off A. In the case at 
issue the relevant poles (the rightmost ones) are found at z = 2, z = 1 and z = — 1. The first one is a simple 
pole due to the divergence of Cm at z = 2. This induces the strongest divergence in E reg (the one that is related 
to the volume (surface) and which in our two dimensional case goes as A 3 ). The next divergence arises from the 
pole of Cm at z = 1 and goes as A 2 . The last divergence stems from a pole at z = —1. This pole is in general a 
double one because both M [g, z + 1] and Cm have poles at that point. This fact means that by properly applying 
the Cauchy theorem, the integrand at z = — 1 determines the finite part of E reg , and also a divergent piece which 
goes as the logarithm of A. Let us give a more detailed account of these remarks. From the hypothesis we have 
stated, we may write 

U[g,z] = - + Fin [M [g, z = 0]}+ 0(z) 
CW = + 0((*-2)°) 

CW = ^^ + o((z-if) 

= R^[C '1] (L7) 

z + 1 

where the symbol Fin means the extraction of the finite coefficient from the Laurent expansion of a function. 



Now we may use these expressions in (1.6) to get that, apart from terms which go to zero for large values of the 
cut-off A, the regularized expression for the vacuum energy is 

E reg (A) = X - (A 3 M [g, 3] Res [C, 2] + A 2 M [g, 2] Res [£ 1] 

+ In ARcs [C, -1] + Fin [M [g] , z = 0] Res [C, -1] + Fin [0 z = -1]) . (1.8) 

The conclusion of this analysis is that to prove that the divergent terms in E reg (A) are independent of the 
magnetic flux, it suffices to show that the residues of Cm at z = 2, z = 1 and z — — 1 do not depend on this 
physical parameter. Let us put it another way, if we label the different heat-kernel coefficients Bq, B±, B\, B^ . . ., 
this independence boils down to saying that Bq, B\ and Bz do not depend on the magnetic flux. As for the 



finite part of expression (1.4) and the quantity Fin [C, z — —1], their relationship is quite direct, as it appears 



explicitly in expression (1. 



We shall systematically throw away such divergences which do not depend on a. These divergences would be 
relevant in a study about the bag dynamics (see |ll|), that is, concerned with situations where the bag walls 
are liable to deformation. In our case, we take the bag to be a perfectly rigid object. 

In the free case (i.e. without flux) the eigenfrequencies under our conditions are zeros of J„ Bessel functions 
with integer indices v coming from angular momentum. The solutions for nonzero a have been found in p, M, 
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and basically correspond to an index shift with respect to the free case \l\ — > \l — a\. Since in both cases the 
eigenmodes are zeros of the same type of functions, we shall introduce the following 'partial-wave' zeta functions 
for fixed values of v: 

oo 

Us) = YjZ> forRes>l, (1.9) 



n=l 



where j„.„ denotes the nth nonvanishing zero of J v (see also [|6| gj; discrete versions of the Bessel problem, 
their solutions and associated zeta functions have also been studied in ) . 

When considering the whole problem in a D-dimensional space, one must take into account the degeneracy 
d(D, I) of each angular mode in D dimensions. Therefore, we define the 'complete' spherical zeta function 

oo oo oc 

Cm(«) = « S E d ^ D ^H^D,l),n = ^ E d(A0C(D,!)W. (i.io) 

l=l min n=l i=imin 

l m i n is the minimum value (if any) of I. In the free case v(D, I) = I + D/2 — 1 and the general form of d(D, I) 
(see e.g. |$5|) is d(D, I) = (21 + D — 2) ^jr^ — ^yp" i but this will change when a flux is present. 

Following the programme we have just put forward, the partial wave zeta function for scalars is obtained in 
sect. 2. From this starting point, we construct the complete zeta functions for D — 2 and D = 3 complex scalar 
fields in sects. 3 and 4, respectively, finding their analytic continuations to s = —I. Numerical results for the 
zero-point energy are then discussed. Afterwards, in sect. 5 we study the Dirac field in D = 2, and sect. 6 is 
devoted to the conclusions. 



2 'Partial- wave' zeta function 

Computing the Casimir energy through the calculation of the complete zeta function requires the knowledge of the 



Bessel zeta functions (1.9) at s = —1, while the complex domain where (1.9) holds is bounded by Re s = 1. This 
is a serious difficulty, but we know that Cv( s ) admits an analytic continuation to other values of s. Moreover, in 
refs. |Q and [||| we showed how to obtain an integral representation of this continuation valid for — 1 < Rc s < 0, 
which reads 



C(s) = -sm — 

7T Z 



dxx 



In 



'2-Kxe x l v (x) 



Whenever c^Owe can work out (2.1) by the method explained in 
1 



for -1 < Re s < 0. 
3(J (sec also fill and 



(2.1) 
|), arriving at 



Us) 



4 1 



-v -sin — 



2s V 2 

.8-1 



B 



+ 2-*B| ill,-* ] } „ 



(2.2) 



S N (s, v) + \ P B ( i±I - S ~) i + J^l + J2 J n (s)± 



with 



1, 



1, 



(2.3) 
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In addition 



C{v, x) 



[>>■>■, - : |ln[£(,^)]-g^»|, 

V2~"^(l + x^Y^e-^I^vx), n(x) = v 7 ! + x 2 + In 



l + x/l - r~ 



W a (t) = 

us) = 

U 4 (t) = 



t 5t 6 

8 ~ "24"' 

3t 4 5t 6 



25 i 3 

"384" 
13 1 4 

128 



531 1 5 ^ 221 i 7 



64Q 
71 r 

"32" 



128 
531 t g 

64 



1105 1 9 

" 1152 ' 
339 i To 



32 



565 1 
128 
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(2.4) 



(2.5) 



the key point being that, this way, S N (s, v) is a finite integral at s = — 1. Further, 

J As) = 



JM = 



1b s+3 s 

48 V 2 ' 2 

00 



dx x 



-s-l 



U n (t(x)), t(x) = 



1 



VT 



1 / s + m s 

2 V 2 '~2 



Thus, the expressions for the ^(s^s are easily obtained from the £/„(t)'s in (|2 

the result of the x-integration is like replacing 

U n {t) -> J„( S ) 

1 / s + m s 
^ ^ 2 B — "a 



— ' 

. In fact, since 



(2.6) 



(2.7) 



(2- 



Expression ( [2.2] ) is not valid for v = 0, since it was obtained from a rescaling x — ► vx and subsequent 
application of uniform asymptotic expansions in vx. Moreover, numerically speaking it is little convenient if v 



is very small. An alternative representation valid in these conditions is needed. Starting from (2.1), we subtract 
and add the asymptotic behaviour of the integrand, which gives rise to a logarithmic divergence on integration. 
When doing so, we shall write the large- x expansion of ln[. . .] as follows: 

4^ 2 -l „ / 1 



In 



'2-Kxe x I v {x) 



4^-1 m ( 1 



8Vx 2 + 1 



O 



1 



Thus, the piece we separate can be integrated with the help of (2.7) (m = 1 case) and we are left with 



COO = 

n v {s) = 



S . ITS 

— sin — 

7T 2 



dx : 



n v (s)- 



Av z - 1 



-B 



s + 1 



16 V 2 
^2~nx e~ x I v (x) 



s 
2 

Av 2 



1 



8VX 2 



1 



(2.9) 



(2.10) 



Since the above integral is now finite at s = — 1 we can Laurent-expand without problems around s = — 1, arriving 
at 



Us) = 



l-Av 2 1 



+ ^^(-1 + In 2) + -KA-l) + 0{s + 1) 

an s + 1 87T 7r 



In particular, for v = 0, 7£q( — 1) = —0.00723 and 

Co(s) ' ' 



0.01451 + C(s + 1). 



8tt s + 1 

Note also the vanishing of the s = — 1 pole when v = ±1/2 already explained in ref. 



(2.11) 
(2.12) 
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3 D = 2 bosons 



3.1 'Complete' zeta function 

Next, we go on to the two-dimensional problem. Following ref.[^), one realizes that the eigenmode sum for this 
case gives rise to the following complete spectral zeta function 



Cm{s\u) = a s C\i- a \(s), 



(3.1) 



Since this function has the properties 

Cx(s;a + fc) = Cyw(s;a), k G Z, 
Ow(s;-a) = Cm («;«)) 
it is enough to study it for < a < 1/2. Introducing 

oo 



1=0 



we can write 



(3.2) 



(3.3) 



(3.4) 



Cm{s;o) = Cm{s;o) + Cm(s; 1 - a) 

= aS C|Q|(s) + Cm(s; 1 + a) + Cm(s; 1 - a). 

oo 

Next we insert expression ( |2.2| ) into (3.3) and, realizing that ^^(/ + /3)~ s = Ch(s,P), where Ch stands for the 
Hurwitz zeta function, wc find 



1=0 



Cm(s;P) = j<7 1 a s Cff(s,^) 



S . 7TS 

+a — sin — 

7T 2 



2s 



-D 



s + 1 s 
2 '~2 



2 S ~ 1 J B 



-2-^' 5 + 3 



8+1 



a) }Ch(8-1,0) 



(3.5) 



1=0 



s + 1 s 
2 '~2 



JV 



+J 1 (sX«(s + l./3) + I^^n( s Kif(s + ^/3) 



with the values of ci, 02 and p in (2.3). Taking AT = 4 and Laurent-expanding near s = —1, this may be written 



(m(s;/3) = - 
a 



- 2 Ch(-1,/3) 



35 



Ch(3,/?) 



■£s 4 (-l,i + )3)(I + /3) 



2=0 



— -^(-2,^ + ^(0,^) 
7T In 2 In 2 tt^(/3) 



s + 1 



64 + 16 ^8 



256 



In a — 1 



- U + -ln2 Cff(-2,/3) 



(3.6) 



ic^(-2,/3) + iCff(0^)K^'^ ' » 



6 



Concerning the pole at s = —1 of the complete zeta function, by (3.4), ( [2.11 ) and (3.6), and noticing that 



6f(-2,l + e0+6f(-2,l 



-of. one arrives at 



(M(s;a) = - 
a 



1 1 



+ 0((s+l) ) 



(3.7) 



128 s + 1 

i.e. the residue is independent of a. The reader may check by using the method explained in ref. |24| that this 
independence applies not only to B3 , but also to Bq and Bi . As we have explained in the introduction, this 
property allows us to state that in cut-off regularization the dependence of the vacuum energy on the magnetic 
flux does not appear in the divergent terms. The dependence of the vacuum energy on the magnetic flux is 
completely contained in the finite part of (,m- 

Since we plan to use the same three formulas for calculating the finite parts, it will be necessary to obtain 
Cjj(0,/3) and ( ! H (—2,(3) around (3 = 1. The first is known (see e.g. pi| ) and amounts to 

1 



CMM)=lnr(/3)--ln(27r), 



(3. 



while the second is calculated by numerical evaluation of Ch(~ n > P) fr° m a n integral representation of the deriva- 
tive of £h valid for negative first arguments. 

3.2 Numerical results 



We start by the I = partial wave zeta-functions obtained from ( 2.1 if ) . Since we are supposing a > 0, the results 
will be denoted by 



s + 1 



+ lna + p a +0(s+ 1) 



l-ia 2 



(3.9) 



where the finite parts p a are listed in the second column of table [I] The pole absence for a = 1/2 may be 



a 


Pa 


Pl+a 


Pl-a 


q a 





-0.01451 


+0.01174 


+0.01174 


+0.00899 


0.1 


-0.05971 


+0.04062 


-0.01172 


-0.03081 


0.2 


-0.10771 


+0.07491 


-0.02987 


-0.06266 


0.3 


-0.15778 


+0.11462 


-0.04285 


-0.08601 


0.4 


-0.20932 


+0.15968 


-0.05095 


-0.10060 


1 
2 


_JL = -0.26180 


+0.21001 


-0.05471 


-0.10650 



Table 1: Finite parts at s — — 1 of the involved zeta funcions (for a — 1). Column 2: I = partial wave zeta function 
Ca(s) Columns 3 and 4: Cm(s;/3) with j3 = 1 + a and /3 = 1 — a. Column 5: finite part of the complete zeta function 
Cm(s;o). 



regarded as a consequence of the fact that J\/2{x) oc sinx, and therefore Ci/aO^) — R~ s Cr{ s ) (Cr meaning the 
Riemann zeta function), which is finite at s = — 1 because £#(— 1) = —1/12. Next, we find Cm( s ',P) from (3.6) 
for the corresponding j3 = 1 ± a's. We shall employ the notation 



Cm0;/3) = 



r 



s + 1 



+ In a 



■P/3 



+ 0(s + l). 



(3.10) 
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According to p.q) 



ra = — 

7T 



(3.11) 



(note that in terms of Bernoulli polynomials n,x) = —^^B n+ i{x)). As for pp, we list some of its values 
in columns 3 and 4 of table [p. Using now (3.4) and the above results we get 



Cm(s;«) 



1 



1 



1 



In a 



Qa 



+ 0(8 + 1), 



(3.12) 
1 

~128' 



128 \s + 1 

The already remarked a-indepedence of the resdiue is exhibited by the fact that r a + f\ +a + fi— c 
Values of q a = p a + p\+ a +pi-a for different a's between and 1/2 are given in the fifth column of table |l| (see 
also Fig. 1). Now it would be incorrect to say that the dependence of E reg (A) on a is exactly given by j^Qa- 
We should take into account a factor 2 which stems from the complex nature of the scalar field. In other words, 
the dependence of E reg (A) on a is given by — q a , apart from terms which vanish when A goes to infinity. 



4 D = 3 bosons 

Eq. ([O]) is again considered, but now in D = 3 and with a magnetic flux line diametrically threading a sphere 
of radius a. We make such a gauge choice that the vector potential in spherical coordinates reads 



eA(r) 



r smf 



(4.1) 



The spectrum and eigenfuctions for the associated quantum-mechanical problem have been written down by 
the authors in ref . [Eg . The proof that one must impose regularity at the origin may be carried out as described 
in ref. |lC| . After studying their associated degeneracies, we are able to write the complete zeta function as follows 



oo oo 



Cx(s;a) = a s ^2 ^2 C|m-o|+p+l/2(s) 
p—0 m— — oo 



again, it is apparent that 

for any integer k, and 

It is now an immediate result that 



CM(s;a + k) = C M {s;a) 

Cm{s; l - a) = Cm(s;cx). 
Cm(s;ol) = Cm(s; -a) 



(4.2) 

(4.3) 
(4.4) 
(4.5) 



From (4.3) and (4.4) we may also restrict our study to the domain < a < h, this is a property which we proceed 



to take advantge of in the sequel. Under this restriction we may give an alternative representaion for (4.2): 

oo 

( M (s;a) = a» ]T |J|C|i-«+i/a|(«)- (4-6) 

In terms of 



l=— OG 



(4.7) 



and of the Cm function defined in (|3.3|), Cm{s\ ol) reads 



Cm(s; a) = Cm s; tt + " + Cm s; ~ + a + Cm s; 



1 



(4.8) 



Next, let's consider the relation between them and the new zeta function 

oo 

C^(s;/3)^a s ^(/ + /3)C ;+/3 ( S ). (4.9) 

1=0 

This has the advantage that Cm{ s > P) can be immediately found from known material. The case without magnetic 
flux has been already studied for D = 3 in (2^, [50|. Since d(3, 1) = 21 + 1 = 2v(3, 1), formula ( 1.1C ) is now rewritten 

oo 

as Cm (s; a = 0) = 2a s i/(3, i)Cw(3,i) ( s )- Therefore, the expression for Cvj (s; /?) is the one for the Cm( s ) m those 
works but for the simple replacement 



Cm(s;P) = -CM(s;a = 0)Ml) = (1 + 1/2) — > (/ + /?)} 



(4.10) 



Thus, for TV = 4 subtractions we find 



Cm(s;/3) = - 

a 



35 



Cff(2,/9) 



1 

H — 

7T 



229 1 1 
._--C*(-3,0 + g fcr(-l,/3) 



1 



s + 1 



In a — 1 



j2s A (-ij+m+p? 

229 



1=0 

293 



(4.11) 



:(ln2-VG8)) 



24192 40320 v 
+ (-1 - ^ ln2j Cff(-3,/3) ~ §Cir(-3,/3) 

+ iln2CK(-l,/3) + ^(-l,^)}+0(s + l) 



As a result of the previous definitions, 



Cm(s;/?) = Cm(s;/3)-/3Cm(s;/?)- 



Hence, the complete zeta function (4.S) is conveniently put in the way 



vvil^: n 1 = ( - - a 



1 



Cm [ s; - + a ) - Cm [ s; - - a 



1 



1 



+ Cm s; 7: + a + (m h; 7 - a , 



1 



(4.12) 



(4.13) 



and the necessary knowledge about the objects on the r.h.s. is available. We have already found the residue of 
Cm(s;/3) at s = —1, namely 

Res \Cm~(s;/3);s = -l] = -fp, (4.14) 
where fp is the one in (3.11). Similarly, from ( |4.1l| ) 



Res 



Cm(s;P);s 



1 

(ITT 



229 1 1 
^ - 5 Ch(-3,/3) + -Ch(-1^) 



(4.15) 



(At this point, one can check that 



Res 



Cm(*;/J = i/2);« = -i 



i l 



a 3157T 



= Res 



:Cx(s;« = 0);s = -1 
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as should be). Then, (4.13) yields 



Res [Cm(s; a); s = -1] = 



2 1, o 

3l5-6 a(1 ~ a 



1 - 



(4.16) 



We recall that this is valid for < a < h, and one has to make use of (4.3) and (4.4) to extend it to any value 



in particular, when we extend expression ( 4.16|) to any real a, we have a non-analytic function of a. The theory 
cannot be renormalized. The situation is different if we give up the idea of a purely confining enclosure and allow 
the presence of external modes, but satisfying the same b.c. as the internal ones. Parallelling the steps in ref. |}(| 
for the a = case, we construct the a-dependent complete zeta function for these external Dirichlet modes — say 
Q M ext (s; a). With respect to the internal case, we have the following modifications: 



(4.17) 



V 7T 

and the i/-series undergoes a (/-parity change, which brings about the transformations 



a 2 
P 

U n {t) 

Jl(s),J n (s) 



-CT2, 

(-i) n z/«(f), 
-Ji(a),(-i) n J„(«). 



(4.18) 



The external ^-function ensuing from these transformations takes into account an overall subtraction from the 
Minkowsky space (for instance, the residue of the rightmost pole is negative). From the construction of this 
external function, all the terms contributing to the s = — 1 pole — including a piece proportional to 1/3 (s) - 
reverse their sign with respect to their internal counterparts and 



Res [Cxoxt( s ;"); s = -1] = - Res [Cx(s;a);s 



(4.19) 



as a result of which the net zeta function Cm( s ! a ) + Cm cxt( s ! a ) * s finite at s = — 1 regardless of the a value. The 
same cancellation applies for the resudues at s = 1 and s — 3. Such a cancellation is typical of odd D's, and does 
not happen in D = 2 because the residue receives then a contribution from (s), which maintains its sign. To 
be brief, we only have to worry about the residue at s = 2. It is quite immediate that it is a independent. The 
conclusion is that for a 3 — D Klein-Gordon field defined in both the exterior and the interior region, the whole 
dependence of the vacuum energy on the a parameter is contained in the finite part of Cm( s ! a ) + Cxcxt( s ! a )- 
The a-dependences of the residue r a and of the finite part p a of (Mi s 'i a ) at s = — 1 are depicted in Figs. 
2a and 2b for the internal modes only. The residue is simply formula ( 4.16| ), while the finite parts have been 
obtained through numerical evaluation of (4.13) by the methods described in refs. 30 . Fig. 2c shows the 
inclusion of the external modes, and the net dependence on a of the vacuum energy. Though we do not pay too 
much attention to the absolute figures, but only to the dependence on a (in other words, to the derivative of the 
finite part with respect to a), it is worthwhile noting that the value at a = furnishes us with an opportunity 
to verify our results. We have obtained that the value of the graph at a = is ^0.005634... = 2 • ^0.002817... i.e. 
twice the figure found in |p3| for an ordinary free field, as had to be expected. 
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5 D = 2 fermions 



For D = 2 massless Dirac particles under the influence of the same magnetic field as in sect. 1 and 3, the Dirac 
equation reads 

(i p + e = Q, (5.1) 



with jb = 7 M f AI (7 = cr 3 , 7 1 = ia 1 , 7 2 = —id 1 ), A°(r) = and A(r) as in (1.2), (for previous works where 
C-function techniques are studied in fermionic systems see [36|, p7t). 

The boundary conditions that we choose on the bag, given by the circle r = a, are those of the M.I.T. bag 
model 

-i = * 

where n stands for the normal vector. 

It has been remarked in refs. j|, ||] that in this problem it would be too restrictive to impose regularity at the 
origin for the modes. If one imposes regularity the result is that the domain of the operator is not dense and, 
consequently, one loses self-adjointness. Making ^(x, t) — ip(x) e~ tEt , let us note the space-dependent part of a 
particular mode by 

X 2 (r) e lv 

It is easily seen that if one demands regularity for the modes characterized by m = — [a] — 1, then one is left only 
with the trivial solution "J = 0. As we advanced in the introduction, we have followed the analysis which was set 



ip($(r, ip)) 



forth in 1 1 ; the outcome is that for the particular value of m = — [a] — 1, one should choose the solution with 
regular x 1 f° r positive a, and the one with regular \ 2 when a is negative. In any case this amounts to picking an 
element among a family of possible self-adjoint extensions for the Hamiltonian under the b.c. in question. Then, 
the whole set of Hamiltonian eigenfrequencies consists of: 

1. the fc's satisfying f u (ka) = J 2 (ka) — J 2 +1 (ka) = with v = {a} — 1 if a > 0, and with v = —{a} otherwise. 

2. (a) the fc's satisfying // + { Q j(fca) = 0, for / = 0, 1, 2, . . . 
(b) the fc's satisfying /j + i_{ ct }(fca) = 0, for I = 0, 1, 2, . . . 

where {a} denotes the fractional part of a. 
It is then adequate to define 

00 

£(s) = K,n, for Re s > 1, (5.2) 

71=1 

where A„ ; „ means the nth nonvanishing zero of f v (X). Now we may write down the fermionic zeta function as 

&(*) = CiaiW+C2mM, (5-3) 

where 

Cl M (s) (*)+0(-a)Cl {Q} (*), 

OO OO 

£m(») = E Cw+n(*) + £CL W+ „(s). (5-4) 



11 



It follows that the Casimir energy will have to fulfil the equalities 



E c (-a) = E c {a), 
E c (a + k) = E c (a), for a > 0, k = 1, 2, 3, . 
E c (a-k) = E c (a), for a < 0, A; = 1, 2, 3, . 



(5.5) 



(compare with relations (3.2)). Since we have now this sort of periodicity when shifting a by integer values, it 
suffices for our study to take < a < 1. 
With the help of the auxiliary object 



(5.6) 



1=0 



(analogous to (3.3) for bosons) we are able to express the complete zeta function as 

£ M {s;a) = a s Ci_!(s) + C^( S ;a)+^;( S ;l-a)^ 

= « s [Ca-i (*) + C («)] + & (*; 1 + a) + C f M& 1 - «) 



(5.7) 



(For the numerical methods to be applied below, the second form proves to be more suitable around a = 0). 



Starting from the partial wave zeta function (5.2), we make use of the technique described in pq , |27j ] and find 
an analytic continuation to the domain — 1 < Re s < given by the integral representation 



fir \ s . ns 

7T Z 



/>oo 

/ dx^- s_1 ln{TOe- 2x [^(x) , for -1< Rc s < 0. 

Jo 



(5.8) 



For zy ^ 0, and by a subtraction method similar to the one applied in the bosonic case, we obtain the more 
convenient form 

f 1 



_ 8 . TTS 

v — sin — 



1 s + 1 s_ 
s V 2 ' 2 



-S { N {s,v) 



2s 



-B 



s + l s 
2 '~2 



2 S B 



-5 



(5.9) 



s 
2/ 



with 



where 



U{{t) 

ufa) 

Ufa) 
Ufa) 



Slf(s,v) = 
£ { (v,x) = 



In(l-f), 
_t_ f_ 
4 + 12' 

i 8 



1 In [£ f (,,*)]- £ 



n=2 

U l n (t(x)) 1 



n=0 



/2^(l + ^) 1 ^ e -^)]2l [7 2 M+7 2 +iM]) 



(5.10) 



8 , 8 



5i d 

lp + 8" " 
t K 17t 5 

32 + 1~28 



9t 5 



23i 7 



3f h 



179t 9 



(5.11) 



¥ 165t 7 64 37i 8 8 327 t 9 ~ 57 1 10 



576 

7 J-9 



128 



64 



128 



32 



179 1 1 
128 



71t i: 
64 
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J\(s) = ±B 



2 

jfo) = I dxx-^U'Mx)). 



and 



o 

In order to include v = or close values, we also find the alternative representation 

2 



(5.12) 



Cl(s) = 



1 



1 



1 



1 



---li/ 

1 7T 



dx : 



which is like ( 2.11 ), for bosons. Note the vanishing of the s 

C=-i/ 2 ( s ) = W 2 )^ 25 - l)Cfl(s) is finite at s = -1. 

The analogue of expression (|3.6| ) for the fermionic case is 

/ 97 



In 2 
1 pole for t/ 



^i(-i) 

7T 



C(s 
1 



1), 



(5.13) 



1/2, which happens because 



Ic ff (-2,/3) + ^Cif(0,/9) - 



V 20160 



- 

256 / 



fe( 2 ,/?) 



( 13 



+ 



fi£(-l, i + j8)(Z + 0) 



-i + f + ^-&(-2,/3)-&(-l,/3) 



1 

24 



In 2 



In 2 V(/3) 7n/»(/3) 



a + 1 



In a — 1 



35tt 



V 20160 32768 



Or (3,0) 



(5.14) 



24 



1 



128 



- (2 + In 2) Ck(-2, /?)-(! + In 2) /?) 



-C^(-2,/3) - ~ i C^(0,/3) + 1) 



Actually, with this plus (p . 7|) and (5.13) we realize that 



( f M (s:a) = - 
a 



1 

64 



•ln( 



(5.15) 



i.e. after adding up all the contributions, the residue of the resulting pole at s = — 1 is independent of a, like the 
residues at s = 1 and s — 2, and, in consequence, the dependence of E vac (A) on a is given by —-^q a (remember 
the minus sign associated to Dirac particles, see for instance ref. [|l9|), where q f a is the remaining finite part of the 
^-function once the pole has been removed and does — predictably — depend on the magnetic field. Numerical 
values are given in table [| The fact that 

Qi = 9oj comes from the equality Ca^( s ; 0) = Cm( s i 1)j which is in the 
end a consequence of the modified Bessel function identity 7_ n = /„, n = 1, 2, 3, . . .. The values of q f a are shown 
in Fig. 3a. For the sake of clarity we also give a plot for an extended domain of a which illustrates the particular 
periodicities of the fermionic case (5.5), see Fig. 3b. 



6 Ending comments 

A scalar Klein-Gordon field subject to Dirichlct boundary conditions and under the influence of an external 
magnetic field producing a single flux line has been studied in two- and three-dimensional spaces. We have 
obtained a nontrivial effect in the dependence of the vacuum energy on the flux which would be invisible (within 
the order of our approximation) without the presence of a finite-sized bag. Considering only the internal field 
modes in the D — 2 case, we arrive at the conclusion that the vacuum energy undergoes a finite variation when 
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a 







-0.00583 


0.1 


+0.05603 


0.2 


+0.07735 


0.3 


+0.07753 


0.4 


+0.06656 


0.5 


+0.05045 


0.6 


+0.03314 


0.7 


+0.01733 


0.8 


+0.00484 


0.9 


-0.00312 


1 


-0.00583 



Table 2: Finite part q a of the fermionic complete zeta function for varying reduced flux a. 

the magnetic flux is changed. It is interesting to note that around a — 0, the vacuum energy decreases when the 
flux grows. In this sense, the system would seem to energetically favour the presence of such fluxons. 

In D = 3 the flux line is diametrically threading a sphere and we have quite a different resulting picture. 
Taking just the internal modes, we see that it is not true that the divergent pieces are independent of the magnetic 
flux. In fact we have seen that the coefficient giving the logarithmic divergence in A (or if you prefer, the residue 
of the ^-function at z = —1) depends non-analytically on a. This is also the case for the coefficient giving a 
A 2 divergence, associated to the residue at s — 1. The inclusion of external modes dramatically modifies the 
situation. Their associated divergences exactly cancel those from the internal part, except the one arising from 
the pole at s — 2, but this piece does not contain any dependence on a. At a = 0, our result agrees with the one 
found in ref . |23| . Around this point the system seems to oppose to the growth of the kind of fluxons we have 
pictured, in the sense that some amount of energy must be provided. 

Fermions in D = 2 have also been considered. While the bosonic energy was periodic in a with period= 1, 
the fermionic one is — not too surprisingly — periodic with the same period only on each real semiaxis separately 
as shown in Fig. 3b, where q f a is represented. The divergences of E reg (A) are independent of a, with the 
same transparency in the physical interpretation of the result as in the D — 2 bosonic case. For values in a 
neighbourhood of a = 0, the energy is seen to decrease as the absolute value of the flux grows. So we have that 
the fermionic case in D = 2 shares with the scalar one this property. 

Refs. [ p5| include a study of the three-dimensional bag involving gauge (bosonic) and fermionic massless fields 
without external flux. By way of rough comparison with some of the figures obtained in these works we may 
evaluate the ratio between the maximum variation of the vacuum energy for Dirac field and the same quantity 
for a complex Klein-Gordon field. Taking into account that this variation is given by °'°^ 97 for the fermionic 
case, and 0975 for the scalar one, we have that the ratio is 0.41. In other words, the energy of a Klein-Gordon 
field is in this sense more sensitive to flux changes. 
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To finish this work we shall briefly comment how the analysis that we have performed in this article with 
models without a mass term carries over to cases where this term is present. Of course, had we incorporated 
a mass term, the result for the finite contributions to the zeta functions would be different and would call for 
some extra, though feasible, effort (see flT|| ). The divergent pieces would also change, but in a way which is quite 
trivial. In general, the residue of a pole at a point s, would be transformed into itself plus a linear combination 
of the residues of the poles at s + 2k for positive integer fc's, with coefficients given by even powers of the mass. 
In other words, if one finds that in the massless case divergences are a independent, this property carries over to 
the massive case. 

Acknowledgements 

I. Brevik, E.N. Bukina, E. Elizalde, A. A. Kvitsinsky and K.A. Milton are thanked for comments and discus- 
sions. S.L. gratefully acknowledges an FI grant from Generalitat de Catalunya. 

References 

[1] Y. Aharonov and D. Bohm, Phys. Rev. 115 (1956) 485, 123 (1961) 1511, 125 (1962) 2192, 130 (1963) 
1625; 

G. Morandi and E. Menossi, Eur. J. Phys. 5 (1984) 49. 

[2] Ph. de Sousa Gerbert and R. Jackiw, Comm. Math. Phys. 124 (1989) 229. 

[3] Ph. de Sousa Gerbert, Phys. Rev. D 40 (1989) 1346. 

[4] D. Kabat, Nucl. Phys B 453 (1995) 281. 

[5] L.H. Duru, Foundations of Physics 23 (1993) 809-818. 

[6] M.V. Berry, J. Phys. A 19 (1986) 2281; J. Phys. A 20 (1987) 2389. 

[7] M.V. Berry and M. Robnik, J. Phys. A 19 (1986) 649. 

[8] C. Itzykson, P. Moussa and J.M. Luck, J. Phys. A 19, LI 11 (1986); 
R. M. Ziff, J. Phys. A 19 (1986) 3923. 

[9] F. Steiner, Fortschr. Phys. 35 (1987) 87. 

[10] M.G. Alford, J. March-Russell and F. Wilczck, Nucl. Phys. B328 (1989) 140. 

[11] M. Bordag, E. Elizalde, K. Kirstcn and S. Leseduarte, to appear. 

[12] H.B.G. Casimir, Proc. Ron. Ned. Akad. Wetenschap 51 (1948) 793. 

[13] J. Ambj0rn and S. Wolfram, Ann. Phys. 147 (1983) 1. 



15 



[14] A. Salam and J. Strathdee, Nucl. Phys. B90 (1975) 203; 
L.S. Brown and G.J. MacLay, Phys. Rev. 184 (1969) 1272; 
J.S. Dowker and R. Critchley, Phys. Rev. D 13 (1976) 3224; 
S.W. Hawking, Commun. Math. Phys. 55 (1978) 133. 

[15] D.G.C. McKeon and T.N. Sherry, Phys. Rev. Lett. 59 (1987) 532; Phys. Rev. D 35 (1987) 3854; 
A. Rebhan, Phys. Rev. D 39 (1989) 3101. 



[16 
[17 

[is; 

[19 
[20 
[21 
[22 
[23 
[24 

[25 

[26 
[27 
[28 



E.R. Speer, J. Math. Phys. 15 (1974) 1. 

M. De Francia, Phys. Rev. D 50 (1994) 2908; 

M. De Francia, H. Falomir, E. M. Santangelo, Phys. Rev. D 45 (1992) 2129. 
G. Plunien, B. Miiller and W. Greiner, Phys. Rep. 134 (1986) 87. 

W. Greiner, B. Miiller and J. Rafelski, Quantum electrodynamics of strong fields, Springer, 1985. 

S.K. Blau, M. Visser and A. Wipf, Nucl. Phys. B310 (1988) 163. 

G. Cognola, L. Vanzo and S. Zerbini, J. Math. Phys. 33 (1992) 222. 

C.G. Beneventano and E.M. Santangelo, Int. J. Mod. Phys. All (1996) 2871. 

CM. Bender and K.A. Milton, Phys. Rev. D 50 (1994) 6547. 

M. Bordag and K. Kirsten, Heat-kernel coefficients of the Laplace operator on the 3-dimensional ball, 



preprint UB-ECM-PF 95/1, |hcp-th/9501064| ; 

M. Bordag, E. Elizalde and K. Kirsten, J. Math. Phys. 37 (1996) 895. 

K.A. Milton, Phys. Rev. D 22 (1980) 1441; Phys. Rev. D 27 (1983) 439; Ann. Phys. (N.Y.) 150 (1983) 
432. 

E. Elizalde, S. Leseduarte and A. Romeo, J. Phys. A 26 (1993) 2409. 
S. Leseduarte and A. Romeo, J. Phys. A 27 (1994) 2483. 

G.N. Watson, A treatise on the Theory of Bessel Functions, 2nd edition, Cambridge University Press, 
Cambridge (1944); 

N. Kishorc, Proc. Amer. Math. Soc. 14 (1963) 527; 
E.C. Obi, J. Math. Anal. Appl. 52 (1975) 648; 

J. Hawkins, On a zeta function associated with Bessel's equation, PhD Thesis, University of Illinois (1983); 
K.B. Stolarsky, Mathematika 32 (1985) 96. 

[29] A. Romeo, Phys. Rev. D 52 (1995) 7308. 

[30] S. Leseduarte and A. Romeo, Ann. Phys. 250 (1996) 448. 

[31] A.O. Barvinsky, A. Yu. Kamenshchik and LP. Karmazin, Ann. Phys. 219 (1992) 201. 



16 



[32] A. A. Kvitsinsky, J. Phys. A 28 (1995) 1753; J. Math. Anal. Appl. 196 (1995) 947. 

[33] N. Ja. Vilcnkin, Fonctions speciales et theorie de la representation des groups, Dunod, Paris (1969). 

[34] Bateman Manuscript Project, A. Erdelyi et al, Higher Transendental Functions, McGraw-Hill, New York 
(1953). 

[35] P.D. D'Eath and G.V.M. Esposito, Phys. Rev. D43 (1991) 3234. 
[36] K. Kirsten and G. Cognola, Class. Quantum Grav. 13 (1996) 633. 

[37] J. S. Apps, M. Bordag, J. S. Dowker and K. Kirsten, Spectral invariants for the Dirac equation on the D 



ball with various boundary conditons, hep-th/9511060. 



Figure Captions 

Fig. 1. Bosonic zero-point energy Eq in D = 2 for a = 1 (then Ec — q a ) as a function of the reduced flux a. 
Fig. 2. Description of (m(s; ct) — for the D = 3 internal modes only — at s = —1, as a function of the reduced 
flux a: a) residue r a , b) finite part p a , c) comparison of p a for internal and external modes, together with the 
Casimir energy Ec (aEc = Paint +Pacxt), in D — 3, where p a int is the same as in b). 

Fig. 3. a) Finite part of the fermionic zero-point energy Ec — ~^Q a in D = 2 for a = 1 as a function of the 
reduced flux a. b) Same as in a; we have simply enlarged the domain of a considered. 
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Fig. 1 



0.02 



0.00 - 



-0.02 



-0.04 



0.06 



0.08 



-0.10 - 



-0.12 



0.0 0.1 0.2 0.3 0.4 0.5 



Fig. 2a 



0.005 



0.000 



-0.005 



-0.010 



-0.015 




0.0 0.1 0.2 0.3 0.4 0.5 



Fig. 2b 




Fig. 2c 



0.010 



at: 



c net 



0.008 



0.006 



0.004 



Pa int 



0.002 



Pa ext 



0.000 



-0.002 



-0.004 



0.0 0.1 0.2 0.3 0.4 0.5 

a 



Fig. 3a 




Fig. 3b 



0.01 




0.05 1 ' 1 ' 1 ' 1 ' 1 ' 1 ' 1 ' 1 ' 1 

-2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 

a 



